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Abstract 



We use the results of Maher et al. (preceding paper) to construct the matrix 



of j = partial-wave two-body and 2 — ► 3 scattering amplitudes for the scat- 
tering of longitudinally polarized gauge bosons W^, Zl and Higgs bosons H 

> 

£ ■ correct to two loops in the high-energy heavy-Higgs limit ^ » M E » M w 

m 
o . 

CO ■ be completely absorbed into a running quartic Higgs coupling X s = X s (s, Mfj 

and factors which involve small anomalous dimensions and remain near unity 
After diagonalizing the matrix of partial- wave amplitudes, we use an Argand- 



We show explicitly that the energy dependence of the 2 — > 2 amplitudes can 



CD , 

43 ■ diagram analysis to show that the elastic scattering amplitudes are approx- 



imately unitary and weakly interacting for A s 2.3, but that three-loop 



corrections are necessary to restore unitarity for larger values of A s . That 
is, the interactions in the Higgs sector of the standard model are effectively 
strong with respect to the perturbative expansion for X s 2.3. The bound 
X s ^$ 2.3 for a weakly interacting theory translates to a physical Higgs mass 
Mh 380 GeV if the bound is to hold for energies up to a few TeV, or 
Mh < 155 GeV in perturbatively unified theories with mass scales of order 
10 16 GeV. 

PACS number(s): 14.80.Gt, 11.20.Dj, 12.15.Ji 
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I. INTRODUCTION 



In the absence of direct observations of the Higgs boson of the standard model of elec- 
troweak interactions, its mass Mh and quartic coupling A = Mfj/2v 2 remain unknown 
parameters in the theory related through the known value of the vacuum expectation value 
v = {^/2Gf^ ^ = 246 GeV. It is therefore of interest to search for bounds on Mh, and 
to consider the possibility that A is large enough that the theory becomes strongly interact- 
ing in the Higgs sector PH|. The only strict upper bound on M H , M H < 650-800 GeV, 
follows from the so-called triviality bound in theories with elementary scalar fields [|J as 
implemented for the standard model in lattice calculations f|. A detailed analysis |§ of 
two-body scattering in the Higgs sector of the theory shows that strong-interaction effects 
would be essentially invisible for ^/s <C Mh unless Mh is very large, Mh > 4-5 TeV. Such 
a mass would violate the triviality bound, but could be allowed if the standard model is just 
a low-energy remnant of a more complete theory in which, for example, the Higgs boson is 
composite. 

The suppression of the low-energy scattering is a general consequence of chiral symmetry 
and the known constraints on electroweak symmetry breaking In contrast, high-energy 

scattering is not suppressed. The triviality limit on A, A < 3.5, is large enough that the 
standard model can be strongly interacting with the effects visible at energies i/i S> Mh as 
judged by a detailed analysis of the scattering amplitudes in the Higgs sector calculated to 
one loop [HJ§. That analysis indicated that strong coupling sets in substantially below the 
tree-level bound A ~ 8n/3 = 8.38 (or Mh = 1.007 TeV) derived some time ago by Dicus and 



Mathur |10j and Lee, Quigg, and Thacker JTJ] . The question naturally arises as to whether 
the bound would be further strengthened — or substantially weakened — by extending the 
analysis to two loops. The question is especially pertinent because the analysis in || and 
used renormalizat ion-group- improved scattering amplitudes expressed as series in the 
running coupling A s . An effectively strongly interacting theory appears for \ s > 2-2.5 (see 
also inl-O] for less restrictive analyses). The limits are less obvious in the expanded form 
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of the amplitudes restricted to terms of orders A and A 2 with energy-dependent coefficients. 

Our objective here is to extend the previous analyses to two loops using the results for 
the two-body scattering amplitudes for longitudinally polarized gauge bosons W^, Zl and 
Higgs bosons H obtained by Maher et al. |TB| in the preceding paper. We explore the range 
of validity of the two-loop results with increasing A by making an Argand-diagram analysis 
of the diagonalized j = partial-wave scattering amplitudes for two-body scattering in the 
neutral channels W^W£ , ZlZl-, HH, and Z^H. These amplitudes must lie on or inside the 
usual unitarity circle, and we may use major departures from the circle as evidence of a 
failure of the series to converge at order A 3 . Quantitative measures of convergence indicate 
that the renormalization-group-improved perturbation series converges very slowly, if at all, 
for running couplings A s > 2.3. The standard model becomes in effect "electrostrong" in 
the Higgs sector for larger values of X s , even though the perturbative amplitudes are still 
quite small. 

II. THE SCATTERING AMPLITUDES 
A. Background 

It is well known OJ that scattering amplitudes for processes involving longitudinally po- 
larized gauge bosons W L , Zl are enhanced in the limit Mh 3> Mw by powers of M H /M W 
relative to the corresponding amplitudes for transversely polarized gauge bosons. It is there- 
fore natural to study the scattering of Wl,Zl and Higgs bosons when seeking an upper 
bound on Mh or, equivalently, on the quartic Higgs coupling A = M H /2v 2 . The necessary 
calculations can be greatly simplified in the limit of interest, Mh 3> Mw, by use of the 
Goldstone boson equivalence theorem [|H,[3l,|i9l-|2"2|j . This theorem states that the scattering 
amplitudes for processes which involve n longitudinally polarized gauge bosons and any 
number of other external particles, are related to the corresponding scattering amplitudes 
for the scalar bosons w ± , z to which W^, Zl reduce in the limit of vanishing gauge couplings 
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g, g' by 



T(W£, Z L , H,. . .) — {iC) n T{w ± , z,H,.. 



(1) 



The corrections are of orders M\y/y/s and g 2 ,g' 2 . In the renormalization scheme used in 

jg,c = i + o(<?v 2 ). 

Maher et al. |18| used the equivalence theorem to calculate the complete matrix of two- 
body scattering amplitudes for the neutral channels w + w~ , zz, HH, zH in the high-energy, 
heavy-Higgs limit 3> Mh 3> Mw- (The charged channels provide no extra information, 
as we will see.) The high-energy limit y/s ^> Mh, while not of immediate experimental 
interest, introduces the further simplification that only the dimension-four quartic couplings 
contribute to the two-body scattering graphs to leading order in s. The dimension-three 
couplings contribute to the final results only through the renormalization constants, which 
contain contributions from the low-energy region ^/s ^$ Mh- 



Maher et al. [IS] give the renormalized matrix M. of Feynmann transition amplitudes as 



M = ZF R Z 



(2) 



where Tr is a matrix of finite, partially renormalized amplitudes Ar(s, t, u), and Z is a finite 
diagonal matrix of ratios of renormalization constants. With the two-body channels taken 
in the order w + w~ , zz, HH, and zH, Tr has the form 



R 



A R (s)+A R {t) A R (s) A r (s) 

A R (s) A R (s) + A R (t) + A R {u) A R {s) 

A R {s) A r {s) A r {s) + A R (t) + A R {u) 

A R (t) 



\ 



while 



(3) 



Z = diag (l, 1, Z H /Z W , {Z h /Z w ) 1/2 ) . 



(4) 



The renormalization constants Z w and Zh are given in the Appendix in Eqs. (A4) and (A5). 
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We have indicated only the first variable in Ar(s, t,u) in Eq. (3) since this function is 
symmetric under an interchange of the last two variables. Specifically ||18|| , 

A 2 

A R (s,t,u) = -2A 



16tt 2 
( 



16tt 2 ) 2 



16 ln(-i) - 4 ln(-t) - 4 ln(-w) + 2 + 6\/37r 

-192 ln 2 (-s) + 176 ln(-s) + 96^ ln(-s) ^ 
-481n 2 (-£) + 801n(-t) + 24^ ln(-t) 



-48 In 2 ^ 
+60\/51n 



-it 



801n(-u) + 24 v / 3vrln(- 



456^0 + 138^ 



7T 



+240 In 2 4±± - ^ In 2 - 180tt 2 



V 



/ 



(5) 



e _Mr s, 



-72K 2 - C(2) + 180C(3) 

+ 24K 3 + %]f - 162C(2) In 2 

where s = s/Mjj, t = t/Mfj, and u = u/Mjj. The phases are defined so that — s 
while the variables — i and —u are real and positive in the physical scattering region. £(n) 
is the Riemann zeta function, C = 1.01494... is the value of the Clausen function at 
argument 7r/3, and the K's are certain constants evaluated by numerical integration [IE 
K 2 = -0.86518, K 3 = -0.10666, K 5 = 0.92363. 



B. The running coupling and anomalous dimensions 

A renormalization group analysis of the scattering amplitudes indicates that their entire 
energy dependence can be subsumed in the limit of interest (s ^> Mjy) into a running 
coupling As = X s (s,Mh) and factors involving the anomalous dimensions 7™,7z and 7# 
associated with the w ,z and Higgs bosons. The SO(3) symmetry of the theory in the 



w , z sector gives 7 Z = ^y w [113]. If all momenta are scaled by a factor o so that s,t,u 



a 2 s, a 2 t, a 2 u, the scaled and original scattering amplitudes are related by |23| 

M ({a Pi } , A, M B ) = TM ({ Pi } , X s (a 2 s, Mb), M h ) T. (6) 
Here T is a diagonal matrix, 
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(?) 



with 



cxp 



A s (l) /3(A) 



dA 



i = w,H, 



(8) 



and A s (cr) = A s (<r 2 s, M#). The evolution of A s is determined by the equation 



lna 



dA 



(9) 



/a,(i) /3(A) 

It was demonstrated explicitly in [8,9] that the entire energy dependence of the scattering 
amplitudes calculated to one loop could be absorbed in a running coupling \ s = X s (s, Mfj). 
To that order, /? (1) (A) = 24A 2 /16vr 2 and the 7's vanish. The running coupling A s was defined 
in those references, following the definition given by Sirlin and Zucchini for the complete 
gauge theory, to incorporate some naturally occurring constants as well as the usual leading- 
logarithmic dependence on s, 

A / , . , x/s 



X s (s,M H ) = A 



n -1 



I 24 In y~ + 25 - 3V37T 

16tt 2 I M f 



(10) 



The constant A is then related to the muon decay constant and the physical mass of the 
Higgs boson by [25] 



A = M 2 H /2v 2 = G^M 2 H /V2. 



(11) 



The calculation of •j w ,'Jh, and (3 at two loops involves some unusual features because of 
our on-mass-shell renormalization conventions, and is considered in Appendix A. The results 



arc: 



A 



16tt 2 



+ 0(A 3 



lH = (150 - 247^) ^ 



A 



P = 24 



16tt 2 



1 - 13- 



16?r 2 / 
A \ 



16tt 2 



+ 0(A 3 ) 
0(A 3 ). 



(12) 
(13) 
(14) 



Integration of Eq. (§) gives the running coupling \ s , 



1 ' (/3 + /?iA)ln ( x + O(A 2 ), (15) 



X,(a) A s (l 



or with a chosen as the ratio of a/s to the energy scale at which A s = A 



A s (s, M H ) = A 



A ^l-^ 241n#+25-3^ 



i 



(16) 



16?r 2 V 16vr 2 ; v M H 
It is interesting to note that the Landau pole in X s apparently disappears for A > 167r 2 /13 ~ 
12 or Mu > 1214 GeV when the second term in the denominator changes sign. However, as 
we will see later, such values of A are far outside the range in which a two-loop calculation 
is reliable. 

With the results above, the factors Tj in Eq. (||) can be written as 

where 7 0i i and j3 are the coefficients of A 2 in Eqs. (|T2])-(|T4D. With the use of the one- loop 
running coupling in Eq. (|10D, this becomes 

r, = ex P f- 70jl [^J^ + oi^ ■ ( 18 ) 

The X s which appears in this equation is to be interpreted as the running coupling in accord 
with the renormalization group analysis. The scale a will be defined as above so that 



fs 

(3 Q In a = 24 In + 25 - 3v^tt, 



Mm (19) 



or 



^#expf 25 " 3 ^Ul-43#- (20) 
Mh \ 24 J M H ^ ' 

With these conventions understood, 

r i = a- 7l(As) . (21) 

It is straightforward to check that all the energy dependence of the scattering amplitudes 
M. can indeed be absorbed in the running coupling and the anomalous-dimension factors 7» 
as indicated in Eq. (B|) f25fl. Specifically, 



M(s,\) = T(s,\ s )M(cos6,\ s (s,M H ))T(s,\ s ) (22) 

where M. depends on s only through the running coupling, and is otherwise a function only 
of the scattering angle. The anomalous dimensions ^ w and 7# are quite small numerically 
for X s in the range of interest, 

j w = -0.000240A 2 , (23) 
lH = 0.000778A 2 . (24) 

It will therefore be a good approximation in parts of the later analysis to replace the diagonal 
matrix T in Eq. (B^) by 1. 



C. Partial-wave 2^2 amplitudes for j = 

Our later analysis of (apparent) violations of unitarity in perturbation theory will be 
based on the known properties of partial-wave scattering amplitudes. The matrix a 2 ^ 2 of 
2^2 partial-wave scattering amplitudes for angular momentum j is related to the matrix 
Ai of Feynman amplitudes for the various channels by |j,26 



a, 



2^2, 



s) = NAj~N (25) 



where Aj is the properly normalized partial- wave projection of .M, 

= 327 (TIT) ' / I 1 rfcose - M ( s ' cos ^) p j( cos0 )- ( 26 ) 

The momentum-dependent prefactor approaches unity for y/s 3> M#. The matrix N incor- 
porates the symmetry factors which must be inserted for each pair of identical particles in 
the initial and/or final state P,|26f . It is given for channel labels w + w~ , zz, HH, zH by 



N = diag(l, lV2, 1/V2, 1). (27) 

corresponding to two-body states \w + w~), -^\zz), -j=\HH), and \zH) normalized over the 
entire solid angle. The matrix a 2 ^ 2 is related to the S matrix by 



(sr 2 - i)/2i. 



(28) 



We will deal only with the j = partial-wave amplitudes. The scattering amplitudes for 
j > are quite small, and do not give useful unitarity constraints. The j = amplitudes are 
easily calculated, but the results are too lengthy to record in detail. However, to illustrate 
their character, we give the result for the diagonal zH channel: 



s) = a 







(*), 



where a^ 1 is a cubic polynomial in X s , 

2A. S A. A 



xzH, 



s) 



16tt 



+ 



16, lfe2 (* + W3 + «) 

2980 + 360C(3) - 356C(2) - 648fT 5 ^ 



+ 



A., 



A. 



16n \ 16n 2 



V 



-1UK 2 - 640 ln2 - 324C(2)ln2 
+48 In 2 ^tl + 6727rV3 - 1 124^0 
+327^ In 3 + 144^ In 4±± 
-162vr 2 - i(2287r + %ir 2 y/$) 



(29) 



(30) 



D. The 2^3 amplitudes and unitarity 



The unitarity of the S matrix provides a nontrivial check on the calculation of the 
scattering amplitudes. The relation S^S = 1 reduces for the partial- wave amplitudes to the 
relation 



Imaj 



2-+2 



,2^2 



)ta 2- >2 + E(a 2_y a 2_ ! 



(31) 



n>2 



where the generalized sum in the last term includes an integration over the n-particle phase 
space. Expanding the scattering amplitudes in power series in A and equating like powers 
to order A 3 gives the matrix relations 



Im a' 



Im a' 



(o) 

3 

(1) 
(2) 



o, 



,(0)t_(0) 



toaf = af%ea(" + ( R ea«)*af + £ (a ? - )* aj 



2-^3 
tree - 



(32) 



The last equation relates the imaginary part of the two-loop amplitude to the one-loop and 
tree level amplitudes for the 2 — > 2 processes, and the tree-level contributions to the inelastic 
2^3 processes. The diagrams for the 2 — > 3 processes are shown in Fig. 1. 

One would normally expect the contributions of the 2-^3 processes to vanish for s — > oo 
since the relevant diagrams involve dimension-three operators. The graphs in Figs, la and 
lb are suppressed as expected by the propagator of the exchanged particle (proportional to 
1/t or l/u) and vanish for s — > oo as Ins when projected to any fixed j. However, the 
two-particle "jet" in the graph in Fig. lc can have a low mass so that the extra propagator is 
nearly on shell. As a result, the square of this graph is not suppressed in the unitarity sum 
in Eq. (|3~TD after the integration over the 3-body phase space. However, interference terms 
between jet graphs which differ by the exchange of a jet particle and the final particle from 
the 4-point vertex are still suppressed. The only finite contributions of the 2 — > 3 processes 
to the unitarity sum are therefore from terms which have the topology of a cut scattering 



Eye graph in the terminology of [Te[] . The three such graphs for zH scattering are shown in 
Fig. 2. We will consider this set as an example. 

The contribution of any of the diagrams in Fig. 2 to the unitarity sum is given up to 
vertex and symmetry factors by 

|Pa| f dP 3 



4-y/s J Air \si2 — m 2 + iej su — m? + ie 



s>M| (167r) 2 7T J(mi+m 2 ) 2 167TSi 2 

Here A is the triangle function, 



ds 12 A ^X^) (33) 

s\2 — m z ) + e z 



A(a, b, c 



a-b-cf- Abe 1/2 , (34) 



and mi, m,2, y / 5i2, and m are, respectively, the masses of the two particles in the jet, the jet 
itself, and the intermediate particle. The kinematic factor \p a \/Ay/s in Eq. (|33D normalizes 
the 2 —* 3 jet amplitudes (which automatically have j = at 0(A 3 )) to the standard 2-^2 
partial-wave amplitudes used above |27] . 
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The calculations for the graphs in Figs. 2a and 2b are straightforward. Summing over 
the independent diagrams and taking account of the symmetry factors for identical particles 
in the final state we obtain the value 



A 



(4vr) 2 

for the graph in Fig. 2a and 



a®(zH) 



a^(zH) 



(35) 



A 



-9 + 2ttV3) \a { o\zH) 



(36) 



(4tt) 2 

for the graph in Fig. 2b, where a^\zH) = — 2A/167T is the tree-level amplitude for zH 
scattering. 

The graph in Fig. 2c is singular for e = (that is, for Imsi2 — > + ) since the decays 
H — > w + w~ and H — > zz are allowed and the intermediate propagator can be on mass shell. 
A direct calculation retaining the e in the denominator in Eq. ( p2|) gives 

3A 



'tiM 2 h 



(4vr)< 



1 



a^(zH) 



(37) 



where e can be identified in the narrow-width approximation for the H decay as MhTh- 
The tree-level decay width is T H = 3\M H /16n, so the first term in Eq. ( P7| ) is just equal to 
the square of the tree-level zH scattering amplitude. This has already been included in the 
calculation in the first term on the right hand side of Eq. ([H|) since we treated the Higgs 
boson as stable, and should be dropped here — the H always decays. (Alternatively, had 
we recognized from the start that the H does not appear as an asymptotic state, we would 
have had no tree-level amplitude for zH scattering; it would appear as here when the H 
resonance is treated as narrow in intermediate states.) 



With this specification, the combined result for Eqs. (|36"D and (0) reduces to 



A 



;-12 + 27rv / 3)|ai 0) (^) 



-Z 



(i) 

H 



af{zH) 



(38) 



where Z$ is the one-loop contribution to Z H . Upon substituting the results in Eqs. fl3"5p and 
(j38|) into Eq. (|32D , we find that the unitarity condition is satisfied in the channel zH — > zH. 
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Similar results hold in other channels for the remaining non-suppressed 2 — > 3 contri- 
butions. In each case the relevant graphs have the toplogy of cut scattering Eye graphs, 
possibly connecting different 2-body states, and give results which are just the negative of 
the first-order renormalization constant for the intermediate particle in the jet, mul- 
tiplied by a product of a tree-level amplitudes. This is not accidental: after the tree-level 
amplitudes are extracted, the integrals which remain are closely related to dispersion re- 
lations for the renormalization constants. In particular, the self-energy functions satisfy 
once-subtracted dispersion relations 

RelV - Rell m 2 = P - V \ ds — } .. . 39 

Dividing by (p 2 — m 2 ), taking the limit p 2 — > m 2 , and using the definition of Z in terms of 
the physical fields and self-energy functions gives 

z -^i v r „J^L. (40) 

7T J(mi+m 2 ) 2 [S — m 1 ) 1 

The integral in Eq. ([!3|) is already in the form of a dispersion relation for m 2 < (mi+mj) 2 
since e can then be set equal to zero. If m 2 > (mi +m2) 2 , the integration contour is trapped 
between poles at m 2 ±ie, but by moving the contour outside the poles, once above and once 
below and averaging, we pick up a term proportional to e~ 1 lmll(m 2 ) and a principal value 
integral. Once the former is deleted, the integrals which remain give the renormalization 
constants. In particular, the product of the 3-point vertex factors in Fig. 2a and the factor 
A/167rsi2 in Eq. (^) is equal to — ImlT^p (S12), where the minus sign arises from the difference 
between the factor (2i\v) 2 which appears in 11^ and the factor |2iAu| 2 which appears in 
the absolute square of the 2 — > 3 amplitude. The sum of the corresponding contributions 
from Figs. 2b and 2c gives — Imil^. This analysis clearly generalizes to higher orders. 

The final result for the 2^3 contributions to the last of Eqs. (^) — with the H decay 
term eliminated — is of the form 

ar(-EZ (1) K } (41) 
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where is a diagonal matrix composed of the sums of the first-order renormalization 

constants for the final two particles in the amplitude a^), e.g., Z w + Zh for the zH channel. 
With this specification, Eqs. ( |3~2"D are satisfied identically as matrix equations, a useful and 
nontrivial check on the calculations. 
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a 2 -* 2 - - 
a 3 2 



III. ANALYSIS AND CONCLUSIONS 

A. Diagonalization of the partial-wave amplitudes 

The general unitarity relation for the S matrix, S^S = 1, reduces in a basis in which the 
matrix a 2 ^ 2 of 2 — > 2 partial-wave amplitudes is diagonal to the condition 

n>2 ^ 

This gives the familiar constraint that the exact two-body elastic scattering amplitudes lie 
on or inside a circle of radius | centered at {0, |) in the complex plane, a constraint which 
we will use to examine the breakdown of low-order perturbation theory for large values of 
A s . 

The determination of the diagonal 2^2 scattering amplitudes is straightforward. The 
matrix M. of Feynman amplitudes defined in Eqs. would be SO (4) symmetric and 

could be diagonalized in an SO (4) basis except for the difference between the renormalization 
constraints Z w and Zh in the matrix Z. However M. and the partial- wave scattering matrices 
a 2^2 re t a j n the exact SO (3) symmetry of the Higgs Lagrangian JTBJ. We will therefore follow 
0] and diagonalize the scattering matrices using states based on the SO (3) decomposition 
of SO(4). The sixteen possible two-body combinations of w ± ,z,H states break up under 
SO (4) as4®4 = 9©6©1. Because of Bose symmetry, only the symmetric representations 9 
and 1 are allowed for even angular momenta j, and only the antisymmetric representation 6 
for odd j. The SO(4) representations decompose under SO(3) as 9 — ► 5 © 3 © 1, 6 — > 3 © 3', 
and 1 —> 1. The states we will need for even j are |§ 
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9,5 = < 



w + w~ 



w^z 

j=(w+w--zz) , 9,3= \ 
w~z 
w~w~ 

(2w + w~ +zz-3HH), 



w+H 

zH 

w~H 



9,1 



1, 1 = —=(2w + w- + zz + HH), 
v 8 



(43) 



while for j odd the relevant states are 



6,3 = < 



w ' z 

w+w~ ! 6, 3' = < 
w~z 



w+H 

zH 

w~H 



(44) 



Because the S0(3) symmetry is exact, the (9, 5), (9, 3), and two (6,3) representations 
give eigenstates of a 2 ^ 2 to any order in A. However, the two identity representations of 
S0(3), (9, 1) and (1, 1), mix through the S0(4)-breaking contributions of the HH channel. 

The eigenvalues of a 2 "* 2 for any j are easily determined by diagonalizing the 4x4 
scattering matrix for the neutral channels. The charged channels add no new information 
as is evident from the decompositions above. We will concentrate here on j = 0, and 
will use basis states normalized over the entire solid angle, specifically the initial states 
w + w~, zzj\[2, HH/y/2, zH, and the SO(3) basis states Xi,i — l>--->4 defined as 



Xi 


= Xi,i = 


—^(2w + w 


+ ZZ + HH), 


X2 


= Xb,i = 


, — (2w + w 
V2i 


~ + zz-3HH), 






1 ( + - 




X3 


— X9,5 = 


v^ (w w ' 


-zz), 


Xi 


= X9,3 = 


zH. 





(45) 



The transformation O x aO with 
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f J_ i JL 

n/2 >/6 U 



o 



1 

_V1 

2 



1 o 



V 



(46) 




1 

splits off diagonal 9, 5 and 9, 3 amplitudes to any order in A, leaving a 2 x 2 matrix to be 
diagonalized in the X13X2 sector. 

We can determine the eigenvalues of the 2x2 matrix rather easily by transforming to 
the states which diagonalize the matrix to order A 2 , that is, at one loop |J, 



Xi =N{ X2 - V3A X 2), 
X , 2 = N( X 2 + V3A X i). 



(47) 



Here A = | (zjp - Z^) and N = (1 + 3A)- 1 / 2 . In this basis, the off-diagonal elements of 
the 2x2 matrix are of minimum order A 3 and do not contribute to the eigenvalues to the 
accuracy of the two-loop calculation. We conclude, therefore, that the diagonal scattering 
amplitudes a« are given to 0(A 2 ) by the diagonal elements of the transformed matrix, 



cii = fO /T T aOO') 



(48) 



Here 



O' 



1-|A 2 



v^A 



\ 



3 A 2 



V 



(49) 



/ 



■a/3 A 

10 

1 

where we have expanded the normalization factor iV to the relevant order. 

The diagonal elements a\ and a 2 obtained by this construction appear initially as sums 
of terms which involve different factors T(s) from the anomalous dimensions. Since these 
factors are close to unity for the values of s we will consider, it is convenient to expand the 
r's and resum the results to obtain overall effective T's for these two scattering eigenstates. 
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The anomalous dimension factors for CI3 and 04 appear just as given below. The final results 
are best presented numerically: 



where 



ai(s) = a 



71 



167T 

( X, 



-6 + 



A., 



167T 2 



'185.6 + 113.1*) 



a 2 (s) = a 



12 



167T 2 

A., 



167T 

/ A, 



a 3 (s) = a 



73 . 



167T 2 
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(50) 



(51) 



The relative weights with which 7^ and 7^ appear in the effective anomalous dimensions 71 
and 72 can be read off from the probabilities with which the normalized states w + w~, zzj v^2, 
and HH/V2 appear in \i an d X2, Eq- (^). The parameter a is defined in Eq. fl2"0|). The 
a\ amplitude arises primarily from the SO (4) singlet state, and is consistently about three 
times larger than the SO (4) nonet amplitudes. We will use these amplitudes in the analysis 
in the following section. 
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B. Argand diagram analysis 



The diagonalized 2^2 scattering amplitudes a\ and from Eq. ( p0|) are plotted in 
Figs. 3a and 3b as functions of A s . The anomalous dimension factors <r 7i have been omitted 
as they are close to unity for < 100 TeV and the values of X s which will be of interest. 
The corresponding curves for 02 and 04 are quite close to that for 03, Fig. 3b, and are not 
shown. 

It is clear from Fig. 3 that a\ and 03 move away from the unitarity circle rather quickly 
as X s is increased. The contributions of the 2 — > 3 processes to the unitarity sum in Eq. (|32|) 
are quite small, 

A s ^^3-1)1^12^8.5 x 1Q -5 X 3 



,71 



E(ar 3 ) + ar = ^ S (52) 

f^|a 3 | 2 ^2.0x 10- 5 A3 

for scattering from the initial states x'i an d X3- The contributions from 2 — > 4 processes are 
similarly small ||, less than 8.0 x 10~ 5 A^ and 3.2 x 10 -5 A^ for the Xi an d X3 channels for 
\/s < 100 TeV. As a result, the exact 2 — ► 2 amplitudes a\ and 03 must lie essentially on the 



unitarity circle for X s not too large p8[ . The deviations of the calculated amplitudes from 
the circle give a measure of the range of X s in which the calculated amplitudes are reliable. 

The convergence — or lack of convergence — of the perturbation series is illustrated for the 
channel x[ ~ > Xi i n the vector diagram for a\ in Fig. 4. In this figure, we show the zero-, 
one-, and two-loop amplitudes for X s = 2.5 as vectors in the Argand diagram. The complete 
two-loop amplitude + + is the sum of these vectors. It must lie essentially on the 
unitarity circle if the perturbative approximation is to be valid. It is immediately evident 
from the figure that the series is not converging well for X s = 2.5: \a^\ is nearly as large 
as |cjW| which is as large as \a^ + a^\. Furthermore, the imaginary part of the amplitude 
becomes negative for A s ^ 2.6. It must be positve in the exact result. 

We can quantify the incipient breakdown of the perturbation series using several tests 
discussed in ||. We will limit our attention here to familiar ratio tests which quantify the 
observations above about the vector diagram in Fig. 4 and its analogs for the other channels. 
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Further less-general but sometimes more restrictive tests are discussed in |3T 



In Figs. 5a and 5b, we show the ratios \af^ /{af^ + a\ )\ and \a\ /a\ | for the four 
amplitudes a^. For values of X s greater than 2.3 to 3.2 (depending on the ratio and channel 
considered), the ratios exceed unity and there is no evidence of convergence of the series 
from either the ratio of the two-loop amplitude to the previous partial sum, or from the ratio 
of successive terms in the series. The strongest limits, X s < 2.3-2.4 for ratios less than unity, 
come from the X3 channel. This channel is associated with the 5 representation of SO (3), 
and does not involve HH scattering. We also note that Im is negative for X s > 2.26. 

We will adopt the value A Sjmax = 2.3 as the maximum value of X s for which the pertur- 
bation series in X s may reasonably be said to converge at two loops for energies y/s < 100 
TeV. At higher energies, the anomalous- dimension factors in Eq. ([5(]) begin to affect the 
magnitudes of the eij. These factors can be treated in two ways. If the expression for M. is 
used in the form in Eq. (|22| ) given by the renormalization group, the anomalous-dimension 
factors T multiply each of and divide out in the ratios considered above, exactly 

for CI3 and 04, and approximately for the average factors in a\ and ci2. In this approach, the 
limits above are unchanged. An alternative procedure is to expand the anomalous dimension 
factors. This adds an energy-dependent term to a^ 2 \ In this case the ratio bounds from a 3 
and 04 are strengthened and those from a± and CL2 are weakened. For example, \a^/a^\ 
now reaches unity for A s = 1.85 for ^/s = 10 16 GeV and Mu = 155 GeV. The change from 
X s = 2.3 becomes irrelevant when Mh is determined by inverting Eq. ([16]) at the high energy 
scale. A general limit A Sjmax = 2.3 therefore seems reasonable. 



C. Conclusions 

The restriction X s < A s max = 2.3 is a quantitative condition on the running coupling 
which must be satisfied if low-order perturbation theory in A s is to give a good description 
of high-energy w ± ,z,H scattering. As illustrated above, the perturbation series for the a« 
show little or no sign of converging at two loops for X s > A Sjmax . This of course does not 
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preclude at least asymptotic convergence for large X s . It is simply that perturbation theory 
breaks down as a useful tool, and the theory becomes in that sense strongly interacting even 
though the scattering amplitudes are quite small. Nonperturbative methods are then needed 
to investigate the problem. However, it is quite useful to bound the region in which low-order 
perturbation theory can be trusted since that is the method used for most phenomenological 
calculations. 

The limiting value A Sjmax = 2.3 for the running coupling in a weakly interacting or 
perturbative theory translates directly into an upper bound on the mass of the Higgs boson 
in such a theory once it is decided at what energy the bound is to be applied. The running 
coupling \ s = X s (s,M H ) is plotted in Fig. 6 as a function of M H for various values of \fs. 
Different choices for the upper bound on the perturbative region correspond to different 
horizontal lines on this plot. The line for A Sjmax = 2.3 is shown. 

Since there is no experimental evidence at the present time for any breakdown of the 
standard model, it is plausible — and is generally assumed — that it will remain valid up to 
at least a few TeV. If it remains valid up to a (minimal) value y/s = 5 TeV and we apply 
the bound X s < 2.3 there, Fig. 6 shows that the Higgs sector of the standard model will 
be approximately unitary and weakly interacting at two loops only for Mh < 380 GeV. 
This is signifcantly below the values Mh ~ 1 TeV used in a number of phenomenological 
studies done using tree-level amplitudes. These results should be reexamined. The limit is 
also below the nonperturbative "triviality" bound Mh ~ 650-800 GeV obtained in |j| . The 
existence of a Higgs boson with a mass in the nonperturbative or strongly interacting sector 
between our bound and the triviality limit is not precluded, and would be quite interesting. 
Calculations would be difficult. 

The bound on Mh in a perturbative theory becomes much stronger, Mh ^ 155 GeV, if 
the standard model is assumed to hold up to a typical unification energy of order 10 16 GeV. 
This result may be altered slightly when other couplings are included in the renormalization 
group equations, but the Higgs boson will still remain light. 

Finally, as noted elsewhere |9|,|26[|, the limit on the range of validity of perturbation theory 
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can become a real upper bound on Mh if the standard model is embedded in a specific unified 
or dynamical model, which must remain perturbative up to some unification or dynamical 
energy. If this is the actual situation, the Higgs boson cannot be much more massive than 
380 GeV in theories such as technicolor with a low mass scale, and will be much lighter in 
typical unified theories. It will then be accessible in future experiments. 
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APPENDIX: CALCULATION OF 7w , 7H , AND AT TWO LOOPS 

Because of our on-mass-shell convention for the renormalization of the scalar theory, we 
will determine (3 and the 7 's to two loops using the Callan-Symanzik definitions for these 
quantities p3| , 



Here M# and A are the renormalized or physical Higgs boson mass and coupling. M and 
Ao are the corresponding bare quantities. Since the Goldstone boson mass m w is fixed at 
zero, Mh is the only mass in the physical theory. 
The bare coupling A is given in JXS[] as 
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to two loops. Here £ = 47r/z 2 /M^, and /x is the mass scale introduced to keep the coupling A 
dimensionless in d = 4 — 2e dimensions. The renormalization constants Z w and Z# are fl8 
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Finally, Mo, which we will not need explicitly, is given in terms of the self-energy functions 
11° and U° H for the bare fields calculated in @ by p| 



M 2 = M^-Ren° H (M|) + n°(0). 



(A6) 



We have calculated (3 and 7's using two different approaches. The first is based on the 
expressions for the Z's and Aq above. These give the first few terms in series of the form 



00 3 



A = /x 2e A 1 + J2 E a i> 



e\3 



j=l k=0 



^ = 1 + HZ) c ^ 



(A7) 
(A8) 
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Mh and A are regarded in Eqs. (|A1|) and (A2) as implicit functions of the independent 



variables M and A p3; /i and e are fixed. Differentiation of the equation for A gives 



dX _ 8X dX £ dX dM 2 H 
" ~ ~dX dMI ~ Mj ~dtf 3M^' 



(A9) 



or from Eq. 



P = lim2f 



dXo/dC 



(A10) 



s dA /<9A' 

where Ao is defined in the last expressions by the series in Eq. (|A7|) . Substituting the series 
into Eq. ( |A10|) , expanding to order A 3 , and taking the limit e — > gives 

(3 = A 2 (2ai 5 i + 4Aa2,i — SAai^a^o) 
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A potentially singular term of the form 

x 2 e e 



(a>2,2 — a i,i 



(A12) 



vanishes identically in this construction because 02,2 = afi- The result is therefore finite 
and independent of the arbitrary mass //, as expected. 
A similar calculation gives 



= A 2 (a 1]1 ci i o — 202,1), 



(A13) 



hence, 
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We note finally that the mass dimension jm, defined as 



(A14) 
(A15) 
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is connected to f3 and ^ w through the relation A = (Mfi / Mf[Z w ) A. This gives 

7m = /5/2A - 7™ 

= 12 i^- 150 (i^) 2 - (A17) 

A second approach to the calculations above uses the observation that the introduction 
of the mass \x in the course of dimensional regularization is unnecessary. With our on-shell 
renormalization scheme, the theory contains a single mass Mh which can be used to scale 
the coupling in the transition to 4 — 2e dimensions. With the choice fi 2 = Mjj/An, the 
quantity £ disappears from the theory, and 

^A(§)'( 1+ gx>4). (A 18 ) 

oo j—1 \j 

2 = 1+EE^T (A19) 

j=l k=0 e 

Using the Callan-Symanzik definitions, we find that 

' = S3« e > s !ai^ (A20) 

and 

1 F)7 

7< = lim -(3(e)^. (A21) 

It is straightforward to show that the limits of these expressions are identical to all orders 
in A to those obtained from Eqs. ( |A10| ) and ( |A13 ). They are also formally identical to the 



expressions obtained in mass-independent subtraction schemes with dimensional regulariza- 



tion [30], and imply in the same way recurrence relations for the coefficients of powers of 



1/e greater than the first p3 
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FIGURES 

FIG. 1. The three distinct topologies for the tree-level 2^3 processes. The exchange graphs 
(a) and (b) are suppressed in the 2^3 cross section or unitarity sum for s — > oo. The "jet" graph 
(c) is not suppressed because of the contributions of low-mass pairs to the jet. 

FIG. 2. Inelastic 2^3 diagrams with the topology of cut Eye graphs for initial zH states. 
The terminology is that of Maher et al. [18]. We show: (a), the cut £2 graph which contributes to 
— Zw\cLzh\ 2 ; (b), (c), the cut £4 and graphs which contribute to — -^^|a^-| 2 - 

FIG. 3. Behavior of the j = eigenamplitudes for , z, H scattering shown as functions of the 
running coupling A s = A s (s,Mjj): (a) a\ and (b) 03. The factors associated with the anomalous 
dimensions are very close to unity for yfs in the TeV range and have been neglected. The curves 
for d2 and 04 are very similar to that for 03 shown in (b). 

FIG. 4. Vector diagram showing the decomposition of the eigenamplitude a\ into tree-level, 
one- loop, and two-loop contributions a^°\a^\ and a( 2 \ The dashed curve is a segment of the 
unitarity circle. Corresponding values of the running couplings A s are shown along the real axis for 
a,(°\ along the upper solid curve for +a^\ and along the lower solid curve for a' ' + a ( - 2 \ 

FIG. 5. (a) The ratio of the magnitude \a^ | of the two-loop contribution to the j = scattering 

amplitude Oj = o' '+oW+o( 2 'h to the magnitude | a^ ^ + a,^ 1 ^ | of the complete amplitude at one 

loop, shown as a function of the running coupling A s . (b) The ratio |a( 2 )|/|a^| of the magnitudes 
of the two- and one- loop contributions to the scattering amplitude as a function of A s . Requiring 
that each term of the perturbation series be smaller than the previous term (or previous partial 
sum), gives a limit A s < 2.3 for the X3 —> X3 channel. 

FIG. 6. The relation between the running coupling A s and the physical mass of the Higgs 
boson Mjj for various values of the energy at which the unitarity constraint could be applied. 
The horizontal line is at A s = 2.3 which we adopt as the maximum value of X s for which the 
perturbation series in A s may reasonably be said to converge. 



28 



